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Abstract 

^ ! We compute the exact partition function of the 2D Ising Model at critical temperature 

0\ [ but with nonzero magnetic field at the boundary. The model describes a renormalization 



group fiow between the free and fixed conformal boundary conditions in the space of 



CN ■ boundary interactions. For this fiow the universal ground state degeneracy g and the full 



boundary state is computed exactly. 



p • 1. Introduction 

(D : 

^ ■ A conformal field theory (CFT) [1] is understood to be the ultraviolet (UV) limit of 

a class of general quantum field theories (QFT). It is associated with fixed manifolds of 
I the renormalization group (RG) fiow in general QFT. And general QFT's are understood 
^ \ as deformations of CFT's and describing RG trajectories of their corresponding unstable 
manifolds. More precisely, they are interpreted as appropriate relevant perturbations of 
appropriate CFT's [2]. For this interpretation to be viable, it should be possible, at least 
in principle, to compute infrared behavior of a QFT from its UV definition. In 2D, for 
certain "integrable perturbations" of CFT's, this is made possible by the fact that the 
scattering (S) matrix is exactly calculable [2]. 

The S matrix is understood to encode all information about its underlying QFT. One 
partially successful way of reconstructing the QFT from the S matrix is the Thermody- 
namic Bethe Ansatz (TBA) approach [3-6]. In particular, the method yields a c function, 
seemingly having all the characteristics of Zamolodchikov's c function [7]. 

All the above understanding has been and is still being extended to QFT with bound- 
ary. CFT on manifolds with boundary is studied in [8-10]. It is shown in [9] that, corre- 
sponding to each conformal family in the (bulk) CFT, there is a boundary condition which 



t email: robin@physics.rutgers.edu 



1 



preserves conformal symmetry. Corresponding to each such conformal boundary condition, 
there is a boundary state containing all information about the boundary condition. This 
boundary state belongs to the Hilbert space of states in the Hamiltonian picture based on 
a spacetime coordinatization in which the boundary is a spacelike curve. In [9] it is shown 
how to classify and compute boundary states. 

The understanding of CFT with boundary is extended to integrable field theory on 
manifolds with boundary in [11]. It is shown that, for a model which is integrable in 
the absence of boundary, when defined on manifolds with boundary, there are certain 
boundary conditions which preserve the integrability of the model. The notion of integrable 
boundary condition is made precise in [11] and a general theory to understand and compute 
exact boundary states and boundary scattering matrices for reflections off boundaries is 
formulated. Boundary scattering matrices for several models have since been computed 
[11-16]. 

Another aspect of boundary QFT and critical phenomena is the universal ground state 
degeneracy factor g [17]. In [9] it is shown how to calculate g for boundary CFT's. In 
[17], g is conjectured to be monotonically decreasing from a less stable to a more stable 
fixed point along the RG flow of boundary interactions. However as of now, there exist 
no fundamental definition of the "gf function" in terms of fundamental quantities of the 
boundary QFT like the correlation functions, and no proof for the "^f theorem" conjecture. 

In [19] and [20], the Thermodynamic Bethe Ansatz (TBA) method has been used 
to compute the boundary contribution to the free energy of boundary integrable models. 
However, because of certain difficulties discussed in these references, it has so far been 
possible to compute only ratios of the g function using the TBA approach. 

One of the simplest conformally invariant boundary models is the Ising model with 
boundary, studied in [8-10]. In [9] the two essentially different conformally invariant bound- 
ary conditions are identified to be (1) free case: boundary spins are free; and (2) fixed case: 
all boundary spins are fixed to -|-1 or —1. 

In this paper, we study more general boundary conditions breaking conformal in- 
variance: the Ising spins interacting with an external magnetic field at the boundary. The 
model is solved on the lattice in [21,22] (see also [23]). In this paper, we demonstrate a short 
cut way to calculate the exact partition function. In addition, we compute the boundary 
state. We compute the partition function on a long cylinder of length L in one Hamiltonian 
representation and transform it to the cross channel representation. This enables us to 
compute the full boundary state exactly. In what follows, we describe the model, elaborate 
on this strategy and explicitly compute the partition function and boundary state. 



2. Ising Model With Boundary Magnetic Field 

We study the Ising model with its boundary spins interacting with a constant external 
magnetic field Hb- This boundary condition interpolates between the free case when 
Hb — and the fixed to ±1 case when Hb — > ±oo. It thus drives the system from the free 
case in the UV down to the fixed case in the IR, which are the fixed points of this boundary 
RG flow [17]. The model is studied in [11,18]. At critical temperature, the continuum limit 
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of the model is described by the action: 



A = 7^ / dxdy[i/jd^ijj + ijjd;,'ip]+ / dtl-^ipijj+l-aal+ih / dt a{t){e^'ip + e^ip){t) (1) 
27r Jx, Jb 47r 2 Jq 

Here, z = x + iy, z — x — iy are complex coordinates on our manifold V; the boundary i3 
is given by a parametric curve B : z = Z{t) , z — Z{t), t being a real parameter; e{t) = 
^Z{t) , e{t) = ^Z{t) are components of the tangent vector (e, e) to the boundary with 
e{t)e{t) = 1; ipjip are free massless Majorana fermion fields; a{t) is auxiliary "boundary" 
fermionic degree of freedom (see [11,18]) with the two point function: 

{a{t) a{t'))free = ^sign{t - t') , (3) 

and d = ^a; h is the boundary coupling constant (appropriately rescaled external field 
Hb) with the dimension of [Length]~^ . The first two terms in (1) describe the confor- 
mally invariant critical Ising model with free boundary conditions. The third boundary 
perturbation term is the boundary spin operator (see [11]): 

(TB(t) =za(t)[e^ + e^](t) (2) 

The boundary spin operator was identified in [9,10]. Thus (1) can be interpreted as the 
free boundary Ising field theory perturbed by the boundary spin operator as- 

It is easily seen from (1) that the fermion fields ip, 'ip satisfy the following boundary 
condition at each boundary: 

{j^+ix)m-ij^-imt) (4) 

where 

A = 47r/i2 (5) 



and 



m = eHt)^Z{t)), m = eHt)tl;{Z{t)) (6) 



In particular, from (6), we see that the free {h = 0) and fixed {h — > ±oo) boundary 
conditions correspond to: 

ip{t) = V'(t), free case, 

■0(t) = — V'(t), fixed case 

At the left {x = 0) and right {x = L) boundaries of the cylinder (see Fig.l) the boundary 
condition (4) takes the forms: 

(-^ + iX)u)tp{—iy) = {-^ — iX)uj^{iy) at left boundary, and (8a) 
dy dy 



{— — h iX)(jjip{L + iy) = {- iX)(jjip{L — iy) at right boundary, (8&) 

ay ay 



where a; = e 4 and a; = e 4 . 



3. Method Outline 

Quantum Field Theory on the space-time of a cylinder of length L (later we will take 
L to be large) and circle length R (see Fig.l), coordinatized by x and y respectively, can be 
viewed in two alternate Hamiltonian pictures: In the first, one considers the (Euclidean) 
time coordinate to run along the length of the cylinder (i.e. x) and the space coordinate 
along the circle (i.e. y). A Hilbert space of states Hr is associated with any time slice 
X = const across the cylinder with fields ip, ip satisfying periodic (for Ramond sector) or 
antiperiodic (for Neveu-Schwarz sector) boundary condition along the circle. Naturally, 
TCr is the same as in the bulk theory since in this picture there is no spatial boundary. 
The partition function is given by: 

Z = {Bright\e-''''''\Bieft) (9) 

where \Bieft) and \Bright) are the boundary states corresponding to the boundary condi- 
tions at the X = and x = L ends of the cylinder respectively, and Hr is the Hamiltonian 
in this picture. 

In the other picture, the space and time coordinates are interchanged (i.e. x becomes 
the space coordinate and y the time coordinate) and we now have a system of particles 
on a line of length L at temperature T = The space of states Ti-LiBieft, Bright) is 
associated with a y = const curve and these states must now satisfy boundary conditions 
Bieft, Bright at X = and x = L cylinder ends respectively. The partition function is 
given by 

Z± = tr((±l)^e-'^^^) (10) 

where the + and — signs correspond to antiperiodic and periodic boundary conditions 
along y, F is the fermion number operator and Hl is the Hamiltonian in this picture. In 
this picture the free fields ip, ^ admit the following decomposition in terms of plane waves: 

^(z) = ^6fe,e^'='^ , ^ = Y^bk,e-"'^' (11) 

ki ki 

where the sum is over the free particle momenta ki which are constrained to satisfy: 

1 + X{ki) = , where X{k) ^ e^'''^^ (^^^] (12) 

\k + iX/ 

where (12) is obtained by imposing (8a) and (8b) on (11). In (11), 6^^ = b^ki (notice 
that solutions of (12) occur in pairs ±ki) and fe^;, bl^ with ki > {ki=0 does not satisfy 
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(12)) are the usual free fermion creation and annihilation operators satisfying standard 

anticommutation relations {bk„,bl^} = and {&fc„,&fc^} = = {4^,4^}. And, 

from (8a), bki = ~*(|~lx)^'=i' i^iplying that bkiS do not represent independent degrees of 
freedom. 

In the next sections we start with the partition function as sum over the states in 
(10) (described by modes in (11), (12)). We want to transform it so that it represents 
the partition function (9) in the cross channel picture. From this we can compute the 
boundary state. 



4. Partition Function 

The Hamiltonian of our model in the finite temperature field theory picture (x is space 
coordinate and y is time coordinate) is: 

HL^EL+Y^ki bl bk, (13) 

ki>Q 

where ki s are the allowed momenta satisfying (12), and El is the ground state energy. 
The partition function (9) is: 

Z^^^-RE, J^(i±e-^'=') (14) 

Here, as before, the ± signs refer to antipcriodic and periodic boundary conditions along 
the cylinder circle satisfied by the fields ip. At thermodynamic equilibrium the free 
energy F± is given by: 

-RF± = -REl + J2 ^og(l ± ^~^''') (15) 

ki>0 

This can be written as: 

-RF± = -REl + ^ [ dk-^ log(l ± e"^^) (16) 

where the integration contour C = C+ + C_ in the complex k plane and the positions of 
poles of the integrand are shown in Fig. 2. In the above, X' = We write the integral 
along contour C_ as 

^ d/c^ log(l ± y^) = ^ dk^\og{l±Yn)- dfc^^^^^log(l ±Yn) (17) 

where Yr = e~^^. Now we change the integration variable A; — > — A; in the second integral 
on the RHS of (17), and noting that 
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we get: 



(19) 



Here, the contour — C_ is shown in Fig. 3; = and — 1. The last term in (19) arises 
as a result of writing log(YR — 1) as analytic continuation of log(l — Yr) and we have used 



From (16), (19) we get: 



/X' ioo 
dfc^^ = log(l + X) ^ =.-log2 (20) 



-RF± = -REl + j dk-^ log(l ± Yr) + 



oo loo 



(21) 





Here C = C+ + -C_ is shown in Fig.3. Now, 



hl = \y. kii^l^i - ^M) = E (^1^1 - 1) (22) 

fei>0 ki>0 

which when compared with (13) implies El = —\ ki. This, as before, can be written 

fei>0 

as: 

El = ^ I dkk^^ (23) 

which yields the following: 

ioo oo 







In (23), the integration contour is the same as in (16). In (24), the two terms in the 
integral in the RHS are the usual bulk and boundary intensive free energies respectively. 
We denote the boundary term as ey. 
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1 f „ 2AA; 





(25) 



2ex goes as Alog(A^)+ nonuniversal terms which we set to 0. Likewise, we set the nonuni- 
versal term in (24) ^ dkk to zero. In (21), we evaluate the integral in the second term 
by parts and residues and get: 



logS+ = 5];iog(l + X(ia;0), ^i = ^^ l^'^ + L for(+) 



l>0 



R 

27rn 



(26) 



logS_ = log(l + X{iu;n)), ujn = n e Z, for(-) 



n>0 



where the Wick rotated contour C(y is shown in Fig. 4. Putting (21), (24), (26) together, 
we get: 



oo 

-RF± = -2Rex + logE± + ^ ^ dklog{l ± Yr) + 



(27) 



In (27), we recognize the third term to be the Casimir energy on the circle. The (+) sign 
corresponds to the Neveu-Schwarz sector and (— ) sign to the Ramond sector. We denote 
these by E^^^: 



oo 



(28) 



The fourth and fifth terms in (27) are to be identified as related with the ground state 
degeneracy factor g{R). 



\ogg±{R)^^ J ciA;-p^log(l±e-'^^) + ^log2 





(29) 



With this, we have: 



Z± = e-^<^\g^{R)e-^^yE^ 
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(30) 



Notice that S-t — > 1 as L 



where 



oo. g± can easily be evaluated: 
a = 2h?R 



(31) 
(32) 

(33) 



4. Boundary States 

Now, we expect that: 

Z± = Zhh ± (34) 

where 

Zhh = (S^|e-^^«|5^) , = (S^|e-^^«|S_^) (35) 

Here is the boundary state corresponding to boundary magnetic field h. In the large 
L limit, the following expressions are valid: 

2Hh = e-'^<\{0\Bh))' + e-'^<~\{a\Br,)f (36) 

Zh-h = e-^<' {mh)f - e-^^H"^ {{a\B^)f (37) 

In (36) and (37), |0) is the ground state of energy E^^^ of the Hilbert space TYfl and it lies 

in the NS sector, while \a) is the lowest energy state of energy i?^ -* in the R sector of Tin: 
|(t) = (T |0), (7 being the spin operator in the bulk theory. Also, in the above, we have used 



{Q\B_h) = {Q\Bh) and {a\B_n) = -{(T\Bh) (38) 

Moreover we have chosen the phase of |0) to be such that {Q\Bh) and {cF\Bh) are real (see 
[17]). This can always be done. From (30), (34), (36), (37), it follows that: 

mn) = ^^+(i?)e-«^^ (39) 

{a\B,) = i=p_(i?)e-^^^ (40) 

Now we construct the boundary states more explicitly. First we consider the picture 
in which x = time and y = space (it is in this picture that the boundary state \Bh) belongs 
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to the Hilbert space Hr) and write down the mode expansions for the i/j and i/j fermion 
fields. In the Neveu-Schwarz and Ramond sectors we have 



Neveu-Schwarz sector : 



'il^{x,y) = ^[a^^ie + <+ie J (41) 



n=0 



oo 



n=0 



where 



/fe„,j=?^, n.Z (43) 

Here (ci„_|_i) and (f^n+i) ^-^^ the creation (annihilation) operators for right 

moving and left moving free Majorana fermions respectively. In the other sector we have: 

Ramond sector : 

oo 

iP{x,y) = [ane-'="(^+*^) + aj,e'="(^+^^)] (44) 

n=0 

oo 

{l,[x,y) = ^ [a^e-'^"^^-^^) + 4e'="("-*2')] (45) 

n=0 

where 

27rn „ 

" = ^ ^ ^ (46) 

and with similar interpretation for the operators ajj, On, a]^, as in the NS case. 

In terms of these mode creation and annihilation operators, the boundary condition (4) 
can be expressed as the following set of constraints on any boundary state 

aABu) = i^^^a\\B^) (47) 

Here i e Z + ^ if the ai s are in the NS sector and Z e Z if the a/s are in the R sector, and 
h = Prom equations (39) and (40) and mode expansions (41) - (46), the boundary 
state takes the following form (following [11]): 
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n=0 



(48) 

± -j=g.{R)e-^'- exp{ 5^/C(/c,) aj, 4}|a) 



where 



lC{ki) = n{iki) = i^ ^ = iL^ , a = 2h^R (49) 

ki + X I + a 

where Tl{k) = is the "massless boundary scattering matrix" for our model [11]; /C 

can also be obtained directly from (47). Putting (31) - (33), (48) and (49) together, we 
finally obtain the boundary state: 



1 cso 1 

2i^a n-a . ] 

r(a+l) ^ n + a " 

^ ^ n=l 

This expression (50) is the main result of this paper. 

Now, we derive from (50) the boundary states corresponding to free and fixed bound- 
ary conditions. These correspond to the UV and IR limits of the boundary RG fiow as 
discussed earlier. 

free case (h — 0) : 

From (50), we readily obtain, putting a = 0, 



n=l n=l 

The above follows since the creation operators anticommute and since the second term in 
(50) is zero in this case. Hence, 



n=l 

since 



oo 

'/,ee)=expi{^ al^, <+i}(|0)-|^)) (52) 



at at \0) ^i\e) (53) 



2 2 



£ being the energy density operator of the bulk critical Ising model. 
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fixed case (h — > ±00 ) : 



In this case, we have from (50): 




n=l 



t 



1} (|o) + k)) ± 



(54) 



± exp-i{^ al 4} 
^4 1 



Both (52) and (54) agree with well known results of Cardy [9] . 



Conclusion 

The expression (50) for \B±h) is the main result of this paper. In the conformal limits, 
it matches with Cardy's results [9]. Our method of computing the boundary state is the 
same in spirit as Cardy's [9,10] in deriving boundary states in boundary CFT - namely, 
making use of a (modular) transformation to express the partition function in one channel 
in terms of states in the cross channel. An essentially similar calculation can be done for 
the off critical Ising model with boundary magnetic field; we will not deal with it here. It 
would be interesting to extend the method to general interacting boundary integrable field 
theories. 
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